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Abstract. We investigate a possible use of direct photon production in association with a heavy 
^ ' quark in pA colHsions at the large hadron collider (LHC) to constrain the nuclear gluon parton 

distribution function (PDF). This process is sensitive to both, the nuclear heavy quark and gluon 
parton distribution functions and is a very promising candidate to help determine the gluon nuclear 
PDF which is still largely untested. 
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&■ INTRODUCTION 

Parton distribution functions (PDFs) are an essential component of any prediction in- 
CS I volving colliding hadrons. The PDFs are non-perturbative objects which have to be 

Q ■ determined from experimental input and link theoretical perturbative QCD predictions 

cn '. to observable phenomena at hadron colliders. In view of their importance, the proton 

^ I PDFs have been a focus of long and dedicated global analyses performed by various 

groups. Over the last decade, global analyses of PDFs in nuclei - or nuclear PDFs 
(nPDFs) - have been performed by several groups: nCTEQ [1, 2], nDS [3], EPS09 [4], 
and HKN07 [5]. In a manner analogous to the proton PDFs, the nPDFs are needed in or- 
der to predict observables in proton-nucleus (pA) and nucleus-nucleus (AA) collisions. 
However, as compared to the proton case, the nuclear parton distribution functions are 
k^ I far less well constrained. Data that can be used in a global analysis are available for 

^ ■ fewer hard processes and also cover a smaller kinematic range. Here, we focus on con- 

straining parton densities in nuclei using the production of a direct photon in association 
with a heavy-quark jet in pA collisions (all details of the analysis presented here can be 
found in [6]). 

DIRECT PHOTONS 

Single direct photons have long been considered an excellent probe of the structure of 
the proton due to their point-like electromagnetic coupling to quarks and due to the 
fact that they escape confinement. Their study can naturally be extended to high-energy 
nuclear collisions where one can use direct photons to investigate the structure of nuclei 
as well. Concentrating on a double inclusive production of a direct photon with a heavy 
quark allows access to different PDF components. Single direct photons couple mostly 
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FIGURE 1. Left: nPDF ratios R^'' = g''^''''{x, Q)/g''{x, Q) (top) , R^^ = cP/^''(x, e)/cP(x, Q) (middle), 
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7?™ = m;:' (jt;,2)/Mi,(;c,e) (bottom) at g = 50 GeV within nCTEQ (solid black line), EPS09 (dashed 
blue line), and HKN07 (dash-dotted red line). The shaded regions correspond to the ;ic-values probed at 
RHIC (x - 10-') and the LHC (x - lO^^). Right: double ratios R'^^^/Rf and R^'^/R^^ using the same 
nPDF sets. 



to valence quarks in the proton or nuclei. By investigating direct photons accompanied 
by heavy quark jets, one gains access to the gluon and the heavy quark PDF. That is 
because, at leading order ^(aa^), the direct photon with a heavy quark arises only 
from gQ — )■ yQ Compton scattering process as opposed to the single photon in which 
case a Compton scattering contribution gq — )• jl competes against a contribution from 
quark annihilation qq — t- yg. At leading order, we see that the initial state for the direct 
photon production with a heavy quark jet depends only on gluon PDF and heavy quark 
PDF where also the latter is often radiatively generated from the gluon leading to even 
stronger dependence of this process on the gluon PDF. 

PROBING NUCLEAR GLUON PDF 



In order to obtain results in hadronic collisions, the partonic NLO cross-section for 
direct photon in association with a heavy quark has to be convoluted with PDFs for 
protons and/or nuclei [7]. For nuclei we show results using the most recent nCTEQ [1], 
EPS09 [4], and HKN07 [5] nuclear PDF sets. Each set of nuclear PDFs is connected 
to a set of proton PDFs to which it reduces in the limit A — > 1 where A is the atomic 
mass number of the nucleus. Therefore we use the various nPDFs together with their 
corresponding proton PDFs in the calculations. 

No current nuclear PDF analysis include a possibility of an intrinsic heavy quark 
component at the input scale and so the heavy quark PDF at a higher scale is generated 
using the gluon PDF where the gluon splits into heavy quark and anti-quark. Therefore, 
the heavy quark PDF follows in shape the gluon PDF rather closely (see Fig. 1) and we 
will concentrate on discussing the gluon PDF in more detail. 

A common feature of all nPDF global analyses is that the nuclear gluon distribution 
is only very weakly constrained in the x-range 0.02 < x < 0.2 from the 2^-dependence 
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FIGURE 2. Nuclear modifications R'^ = gP^'^{x,Q)/g''{x,Q). Left: For lead at Q = 50 GeV. Shown 
are results for nCTEQ decutS (solid, black line), EPS09 (dashed, blue line) + error band, HKN07 (dash- 
dotted, red line) + error band. The box exemplifies the x-region probed at the LHC (y/s = 8.8 TeV). Right: 
nuclear production ratio of 7+ c cross-section at LHC within ALICE PHOS acceptances, using nCTEQ 
decutS (solid black line), nCTEQ decutSgS (dotted black line), nCTEQ decut3g9 (dash-dot-dashed black 
line), EPS09 (dashed blue line) + error band, HKN07 (dash-dotted red line) + error band 



of Structure function ratios in deep-inelastic scattering (DIS) F^"(jc, 2^)/F^(x, 2^), 
measured by the NMC collaboration [8]. In order to compare the various nPDF sets, we 
plot in Fig. 2 the gluon distribution ratio Ri{x, Q) ~ gP'^{x, Q)/ g''{^i Q) as a function 
ofx for a lead nucleus at 2 = 50 GeV. The chosen hard scale Q = 50 GeV is typical for 
direct photon production at the LHC and the box highlights the x-region probed by the 
LHC. 

The fact that the nuclear gluon distribution is poorly constrained is reflected by large 
PDF uncertainty bands of the gluon PDF of the HKN07 and EPS09. Moreover, addi- 
tional uncertainty is connected to the choice of parameterization and other assumptions 
that are an integral part of any PDF analysis. The second source of uncertainties causes 
uncertainty bands of HKN07 and EPS09 to not overlap in some regions and causes also 
the nCTEQ gluon PDF to lie outside of the other uncertainty bands. Also the rather 
narrow and overlapping bands at small x < 0.02 do not reflect any constraints by data, 
but instead are theoretical assumptions imposed on the small-x behavior of the gluon 
distributions. 

The nuclear production ratio 7?^+^' = A ^JXrlipp^ncS) i" a kinematic range 
probed by the ALICE experiment at the LHC is shown in Fig. 2 (right) using several 
nuclear PDFs EPS09 (dashed blue line), HKN07 (dash-dotted red line) and a series of 
nCTEQ fits decut 3 (solid black line), decut 3g3 (dotted black line) and decut 3g9 
(dash-dot-dashed black line) which differ in assumptions on the small x behavior of 
the gluon PDF (details see [6]). For the first two cases the bands represent the nPDF 
uncertainties. Remarkably, there is almost no overlap between the EPS09, the HKN and 
nCTEQ predictions, therefore an appropriate measurement of this process will be able 
to distinguish between the nPDF sets. 

In the ALICE experiment, photons can be identified in the EMCal electromagnetic 
calorimeter, or in the PHOS spectrometer with a somewhat more limited acceptance. 
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FIGURE 3. NLO differential cross-section for 7+ c (left) and y+b (right) production in p-Ph collisions 
at a center-of-mass energy of y^ = 8.8 TeV in PHOS (lower band) and EMCal (upper band) acceptances. 



In Fig. 3, we show results of calculations carried out for p-Ph collisions at the LHC 
nominal energy y/s = 8.8 TeV using the cuts and acceptances of both EMCal and PHOS 
detectors. The differential NLO cross-section is plotted as a function of the photon 
transverse momentum in the y + c (7 + Z?) channel in Fig. 3 left (right) for both PHOS 
(lower band) and EMCal (upper band); the dotted curves indicate the theoretical scale 
uncertainty. 

The total integrated cross-section for the EMCal electromagnetic calorimeter is 
1 19000 pb for 7+ c process and 22700 pb for y+b. 

CONCLUSION 

We have presented a detailed phenomenological study of direct photon production in 
association with a heavy-quark jet in pA collisions at the LHC, at next-to-leading order 
in QCD. The dominant contribution to this process is given by the gQ — )■ 72[+g] 
subprocess. This offers a sensitive mechanism to constrain the heavy-quark and gluon 
distributions in nuclei, whose precise knowledge is necessary in order to predict the rates 
of hard processes in heavy-ion collisions where quark-gluon plasma is expected to be 
formed. 
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